arXiv:1507.05958vl [astro-ph.GA] 21Jul2015 


PROCEEDINGS 

OF SCIENCE 



Multi-wavelength constraints on cosmic-ray leptons 
in the Galaxy 


E. Orlando*", A. W. Strong 6 ,1. V. Moskalenko ', C. Dickinson', S. Digel", T. R. Jaffe^ 
G. Johannesson', J. P. Leahy', T. A. Porter", M. Vidal" 


"Hansen Experimental Physics Laboratory and Kavli Institute for Particle Astrophysics and 
Cosmology, Stanford University, SLAC National Accelerator Laboratory, Stanford, CA 94305, 


U.S.A. E-mail: |eorlando@ Stanford. edu 

b Max-Planck-Institutfiir extraterrestrische Physik, Postfach 1312, D-85741 Garching, Germany 
c Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, 
The University of Manchester, Oxford Road, Manchester, Ml 3 9PL, U.K. 
d CNRS; UPS-OMP; IRAP; Toulouse, France 

e Science Institute, University of Iceland, Dunhaga 5, IS-107 Reykjavik, Iceland 


Cosmic rays (CRs) interact with the gas, the radiation field and the magnetic field in the Milky 
Way, producing diffuse emission from radio to gamma rays. Observations of this diffuse emission 
and comparison with detailed predictions are powerful tools to unveil the CR properties and to 
study CR propagation. We present various GALPROP CR propagation scenarios based on current 
CR measurements. The predicted synchrotron emission is compared to radio surveys, and syn¬ 
chrotron temperature maps from WMAP and Planck , while the predicted interstellar gamma-ray 
emission is compared to Fermi- LAT observations. We show how multi-wavelength observations 
of the Galactic diffuse emission can be used to help constrain the CR lepton spectrum and prop¬ 
agation. Finally we discuss how radio and microwave data could be used in understanding the 
diffuse Galactic gamma-ray emission observed with Fermi-L AT, especially at low energies. 
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1. Cosmic ray leptons and diffuse emission 

Below a few GeV the local interstellar lepton spectrum (LIS) cannot be directly measured, 
because CRs are affected by solar modulation. At higher energies, CR measurements might not 
be representative for the average spectrum in the Galaxy §• During propagation in the Galaxy 
CR leptons are interacting with the interstellar medium (ISM), and are scrambled by the magnetic 
field (B-field). Their spectrum is steepened by energy losses and energy-dependent diffusion, and 
it can also be affected by re-acceleration and production of secondary CRs 0. Hence, tracing 
CR leptons back to their origin and throughout the Galaxy is difficult. However, a way to study 
CR leptons in the Galaxy is to observe the diffuse emission produced by their interaction with 
B-fields via synchrotron emission (in radio and microwave) and with the gas and the interstellar 
radiation field (ISRF) via bremsstrahlung and inverse Compton (IC) emission (in gamma rays). 
These observations, compared with sophisticated propagation models, allow us to gain information 
on the CR propagation in the Galaxy and and their interstellar spectrum. In particular CR leptons 
from ~500 MeV to tens of GeV produce synchrotron emission from tens of MHz to hundreds of 
GHz for a typical B-field of a few /iG 0. The observed synchrotron spectral index depends on 
the spectral indices of CR leptons along the line of sight. Hence this can be used in conjunction 
with direct measurements to construct the full interstellar lepton spectrum from hundreds of MeV 
to TeV. Combining observations of the diffuse emission from radio through gamma rays, with 
the aid of propagation models, then allows us to derive the quantities related to the CRs, such 
as spectrum, density, distribution in the Galaxy, and propagation parameters. This approach has 
certain advantages over studying these data separately by providing independent information on 
the spectra of CR leptons, and hence reducing degeneracy. 

After presenting a summary of our previous studies on CR leptons and diffuse emissions, we 
account here for recent PAMELA CR lepton measurements 0. The derived synchrotron emission 
is calculated and its spectrum is compared with radio surveys and recent observations from the 
9-year WMAP [|5]] and four-year Planck synchrotron temperature maps [jhj]. The gamma-ray diffuse 
emission is calculated and compared to the Fermi Large Area Telescope (LAT) data. 


2. Modeling CRs with GALPROP 


This study uses a numerical model of CR propagation and interactions in the Galaxy, GAL¬ 
PROP 1 . Descriptions of the code can be found in [^. [/]. ||] and references therein. It enables simul¬ 
taneous predictions of observations of CRs, gamma rays and synchrotron radiation. GALPROP 
models were used to analyze and interpret the diffuse gamma rays detected by EGRET, COMP- 
TEL (e.g. §), INTEGRAL (e.g. 0), LAT (e.g. [® 0), and more recently radio 

and microwave observations (e.g. [[|, [0). GALPROP calculates the diffuse emission from pion 
decay, bremsstrahlung and IC, for a user-defined CR source distribution and propagation parame¬ 
ters. Sample gas maps and ISRF are given in JlH]. Calculation of interstellar synchrotron emission 
has been improved [|j] and extended [1_6] to include 3D B-field models, synchrotron polarization 
and a basic model for free-free emission and absorption. Also 3D gas and source distribution mod- 


1 http://galprop.stanford.edu/ 
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Figure 1: Modeled LIS for a plain diffusive propagation model and modulated spectrum compared with 
PAMELA |fl|| and LAT |p7|| measurements, taken mostly during solar minimum. 



Figure 2: Synchrotron spectrum ( black line ) from high latitudes (10° < \b\ < 45°, 0° < / < 360°) for the 
plain diffusive propagation model calculated with the CR leptons as in Fig. [I} Data are: 9-year WMAP 
synchrotron temperature maps ( blue stars) [|]]; Planck synchrotron maps scaled to 23 GHz | [29| | (red point); 
radio surveys described in the text and in [0| (green triangles); 408 MHz from |[2?]| with 3.7 K offset (yellow 
triangle) and reprocessed 408 MHz from ][23|| with 8.9 K offset (green triangle). 


els were recently implemented [JT7] and Johannesson et al. (these Proceedings). For an updated 


overview see [|18|] and Moskalenko et al. (these Proceedings). 


3. Cosmic ray leptons: summary of our previous studies 


In our recent papers [jj. 11, 16, 19], AMS-01 [|20|] measurements at low energies and LAT [21] 


data at higher energies were used to derive the spectral index of the CR lepton spectrum. Unfortu¬ 
nately, CR measurements are not sufficient to constrain the leptons at low energies, since there is a 
strong correlation in the models between the spectral index and the modulation potential. Hence, 


3 













Cosmic rays from multi-wavelength observations 


E. Orlando 


to break this degeneracy we derived the local spectral index at low energies from synchrotron ob¬ 
servations [|j], using a collection of radio surveys and spectral index data over a wide range of 
frequencies and the synchrotron temperature maps from [|^]. This allowed us to constrain the LIS 
at low-energies and propagation models that best fit the observations while avoiding the effects of 
solar modulation. We concluded that the interstellar spectrum turns over rather sharply below a few 
GeV. This is independent from propagation effects, and it reflects the electron injection spectrum 
from the sources. The low-energy falloff in the directly measured electrons, normally attributed 
mainly to modulation, is partially due to the interstellar spectrum. This suggests less solar modu¬ 
lation than usually assumed. We found that it is challenging to describe the observed synchrotron 
spectrum using present diffusive re-acceleration models with a standard Alfven velocity (that works 
fine for other CR species), since the modeled emission from primary and secondary leptons exceeds 
the measured synchrotron below hundreds of MHz. Spatial effects, independently in gamma rays 
with the LAT |[TI|] and in radio and microwaves |[Tt]], were extensively investigated, using different 
CR distributions, propagation halo sizes, B-field and ISM models. We report here on the spectral 
properties only, updating our previous analyses by using recent PAMELA CR measurements, the 
reprocessed 408 MHz map [23], the new 9-year WMAP maps, the recently released Planck maps, 
and gamma-ray observations with the LAT. 


4. Description of our approach 

We use GALROP for modeling CR propagation requiring that the modulated spectrum re¬ 
produces CR measurements. Our method is to extract information on the spectral indexes of the 
injected particles and propagation parameters by simultaneously studying the radio, microwave and 
gamma-ray emissions produced by CR leptons. This allows us to bypass part of the degeneracy of 
the parameters, such as solar modulation, primary electron spectral index and contribution of sec¬ 
ondaries. Specifically, first we derive constraints on the LIS from direct CR measurements. Then 
we calculate the diffuse synchrotron emission in the radio and microwave band and compare with 
synchrotron observations with the aim of tuning the lepton spectrum at low energy to match the 
observations when possible. Then we use this information to generate the leptonic diffuse gamma- 
ray emission. In turn, comparison of the predictions of the diffuse gamma-ray emission with LAT 
observations provide us additional information on CR leptons and propagation. Therefore, gamma- 
ray and synchrotron emission can probe the interstellar spectrum and propagation effects free from 
the effects of modulation. 


5. Observations 


We use the surveys at 22, 45, 150, 408, and 1420 MHz and their zero levels as in [|3j]. Frequen¬ 
cies from 20 MHz to a few GHz account mostly for the synchrotron emission. The emission there 
is produced primarily by electrons with energy below a few GeV. The reprocessed 408 MHz map 
by [23] with better source subtraction is also used. Its offset is taken from [24] as 8.9 K. Thanks 
to the latest release this year p5|], we use here the Planck low frequency separated component 
maps. In addition the seven-year WMAP temperature maps used in [ |T6| | are updated here with the 
9-year WMAP synchrotron maps from 23 to 94 GHz obtained with the Maximum Entropy Method 
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(MEM). We use LAT observations as published in [11], from intermediate latitudes that trace rel¬ 
atively local CRs. The used LAT data include 21 months of observations and Pass 6 DataClean 
event selection. 


6. Propagation models 

We follow 0 and start from the same plain diffusive model that fits synchrotron observations. 
We tune the propagation parameters from B/C, protons, and helium to PAMELA measurements. 
Then we tune the electrons to direct LAT measurements above 100 GeV where the solar modula¬ 
tion is negligible. The modulated low energy spectrum is then tuned to reproduce simultaneously 
PAMELA CR measurements and synchrotron observations, and finally we compare the models to 
gamma-rays to make sure the modeled emission matches data. The local CR lepton spectrum for 
the new propagation model is shown in Lig. |T} Note that PAMELA measures the electron spectrum, 
and its statistical errors increase with energy. At higher energies we reproduce the total lepton spec¬ 
trum as measured by the LAT. The latter combines electrons plus positrons so that the increasing 
positron fraction at high energies is included in our modeling. 


7. Results: synchrotron and gamma-ray spectra 


Ligures [2] and |] show results for the plain diffusion models with the LIS as shown in Lig. [j]. 


The intensity of the B-held has been fit to the 408 MHz map from [23]. The modeled synchrotron 
spectrum for the high-latitude regions reproduces the observations quite well. Below a few GHz 
data from surveys are well reproduced by the model. At a few GHz contamination from free-free in 
this region, not accounted here, is below 25% [[If], p7| |. In microwave we see a good spectral agree¬ 
ment with WMAP maps integrated over the same region of the model, and a good fit to the Planck 
component-separated synchrotron map. This comparison shows that the LIS used here is a good 
representation of the spectrum that produces the synchrotron emission. This suggests that low- 
energy spectral index, solar modulation and the contribution of secondaries are well constrained, 
as previously found by [3|]. The induced gamma-ray emission is shown in Lig. [3] and compared 
to the LAT data for the intermediate latitudes taken from [[TT|] . The model is within the LAT sys¬ 
tematic uncertainties even without the tuning to the data that would account for uncertainties in 
the ISM. Hence, in a first approximation plain diffusive propagation models can reproduce gamma 
rays as well as the diffusive re-acceleration models usually assumed. A detailed fitting procedure 
will be described in a separate work. Lor comparison with published re-acceleration models, Lig. |] 
(left) shows the electron spectrum used in [ ^ | compared with PAMELA measurements. Similar 
re-acceleration models have been used in [ |TT| ] to fit gamma-ray data. However, these models were 
not tuned to fit synchrotron spectral observations. In fact the same figure also shows the calculated 
synchrotron spectrum compared with the same observations reported in Lig. |], and with the B-held 
intensity fit to the 408 MHz map, as above. An excess of the model at the lowest frequencies is 
visible, mostly due to secondary leptons produced by re-acceleration processes. This excess is now 
reduced with respect to [3j] due to a combination of the reprocessed 408 MHz map and the larger 
offset used here. The spectrum does not reproduce the microwave observations as well as the one 
in Lig. [|. since it under-estimates the synchrotron intensity from Planck. 
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Figure 3: Calculated gamma-ray components and spectrum compared with the LAT data from |[TT|l for 
intermediate latitudes (10° < \b\ < 20°, all longitudes), using the plain diffusive propagation model and the 
LIS as in Fig. |TJ Data include statistical ( grey area) and systematic errors ( black bars). Here components are 
not fitted to gamma-ray data, so uncertainties in the gas and ISRF are not accounted for. Spectra for sources 
and isotropic are taken as in [[TTJ], for the most extreme cases reported there. 30% uncertainty is added to 
the isotropic spectrum, following the study in [^] based on various foreground models. The ISRF has been 
scaled to better reproduce gamma rays as found in |[TT|. 



Figure 4: Left, modeled LIS and modulated spectrum for a diffusive re-acceleration model from []28|]. Data 
are as in Fig. |I] The blue line is the model component for secondary leptons. Right: Calculated synchrotron 
spectrum for high latitudes compared with data as in Fig. [| Dashed line is the contribution from secondaries. 


8. Discussion 

In this work we show the feasibility and importance of using multi-wavelength observations, 
especially at radio wavelengths, together with CR measurements, to constrain the interstellar spec¬ 
trum at low-energies, propagation models, and solar modulation effects. We have presented pre¬ 
liminary results and now discuss here the limitations and potential of this approach. 

The exact derivation of the synchrotron maps as obtained by Planck and WMAP have lim¬ 
itations, due to the various assumptions required and degeneracies with separating multiple as- 
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trophysical components including synchrotron, free-free, thermal dust and anomalous microwave 
emissions (AME) ([^, [29||). Disentangling these components in the WMAP and Planck bands re¬ 
quires various assumptions or a priori information such as data-driven templates using ancillary 
data (where each different emission mechanism is at its maximum) or by approximating their spec¬ 
tral properties (see ]|5[ |], [2^. pT|| ). The method used to derive the WMAP MEM temperature maps 
shown here utilized spatial templates of different emission components from external data as priors, 
and the fit was performed to each pixel independently. The spectra for the different components 
were fixed, while the synchrotron spectrum was obtained from polarization observations. The 
resulting synchrotron map is plausibly contaminated by other components, such as AME and free- 
free emission. This could explain the difference with the Planck data ]^9|]. In fact the intensity 
of the Planck synchrotron map is lower than previously found with WMAP , while the AME and 
free-free intensities are higher. This might suggest that the Planck synchrotron intensity could be 
under-estimated, while the WMAP intensity is overestimated. 


The high-quality Planck component maps recently released were combined with complemen¬ 
tary ancillary data: the 9-year WMAP temperature maps and the 408 MHz map, providing an ex¬ 
cellent fit to the data. The synchrotron spectral index from the best model by Jl6[ ] was used as 
the baseline model to allow a better separation of the synchrotron component. Nevertheless, there 
is likely to be a degeneracies between the various low-frequency components, especially between 
AME and synchrotron. Hence, these component-separated maps face uncertainties related to our 
limited knowledge of their spectra and spatial distributions (plus uncertainties in the offset at the 
different frequencies), as discussed in detail in [29]. Note that also the zero levels of the surveys 
are not clearly determined, that could limit our model constraints. [ 24 ] estimated a monopole of 
8.9± 1.3 K in the 408 MHz map [|2h||, which includes any isotropic component (CMB, Galactic and 
extragalactic), which we use for the fit here. However, the determination of the offsets strongly 
depends on the foreground models and data used. In |^]| we adopted a 3.6 K offset, which increased 
the excess at lower frequencies for the diffusive re-acceleration models. Further model-dependent 
studies and forthcoming data at a few GHz up to ~ 15 GHz (e.g. CBASS [32]) will help in sepa¬ 
rating the components and may provide more strict constraints to the lepton spectrum. 

Regarding gamma rays, the sky above 100 MeV is dominated by emission produced by CRs in¬ 
teracting with the gas and ISRF via pion-decay, IC, and bremsstrahlung. Disentangling the different 
components at the LAT energies is challenging and is usually done in a model-dependent approach. 
Uncertainties in the ISM is the major limitation to our modeling and hence in our knowledge of 
CRs, e.g. as found in Jill]. The situation below 100MeV is still unexplored. Extrapolations of 
present models to such low energies predict IC and bremsstrahlung to be the major mechanisms 
of CR-induced emission, which are of leptonic origin. The fact that energies < 100 MeV were 
not deeply investigated after the COMPTEL era, makes it more exiting now with the advent of the 
LAT pass 8 data and its extension to lower energies Jl4| . 33]. However disentangling the different 
components and characterizing the sources below 100 MeV is even more challenging due to the 
relatively large point spread function of the instrument. This highlights the importance of multi¬ 
wavelength observations including the Square Kilometre Array telescope (e.g. [34]) and C-BASS 
] [32| l that will provide supplementary information at radio frequencies. 

Acknowledgment: E.O. acknowledges support via NASA Grant No. NNX13AH72G. GALPROP development 
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